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Recent developments in simulating fundamental quantum field theoretical effects in the 
kinematical context of analogue gravity are reviewed. Specifically, it is argued that a 
curved space-time generalization of the Unruh-Davies effect — the Gibbons-Hawking ef- 
fect in the de Sitter space-time of infiationary cosmological models - can be implemented 
and verified in an ultracold gas of bosonic atoms. 

> 

oo 

O 

>P ' 1. The concept of an effective space-time metric 

^ : .... . . m 

^^ , Curved space-times are familiar from Einstein's theory of gravitationf^ wliere th.e 

metric tensor g , describing distances in a curved space-time with local Lorentz 

"t^ . invariance, is determined by the solution of the Einstein equations. A major prob- 

lem for an experimental investigation of the (kinematical as well as dynamical) 
properties of curved space-times is that generating a significant curvature, equiv- 

' ^ , alent to a (relatively) small curvature radius, is a close to impossible undertaking 

Q ' in manmade laboratories. For example, the effect of the gravitation of the whole 

CJ ' Earth is to cause a deviation from flat space-time on this planet's surface of only 

the order of 10^^ (the ratio of Schwarzschild and Earth radii). The fact that proper 

ij ■ gravitational effects are intrinsically small is basically due to the smallness of New- 

''^ . ton's gravitational constant G = 6.67 x 10~^^ m^kg~^sec~^. Various fundamental 

C^ ' classical and quantum effects in strong gravitational fields, a few of which we dis- 

cuss below, are thus inaccessible for Earth-based experiments. The realm of strong 
gravitational fields (or, equivalently, rapidly accelerating a reference frame to sim- 
ulate gravity according to the equivalence principle), is therefore difficult to reach. 
However, Earth-based experiments are desirable, because they have the obvious 
advantage that they can be prepared and, in particular, repeated under possibly 
different conditions at will. 

The formalism to be described in what follows is aimed at the realization of effec- 
tive curved space-time geometries in perfect fluids, which can indeed be prepared on 
Earth, and which mimic the effects of gravity inasmuch the kinematical properties 
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of curved space-times are concerned. Among such perfect fluids are Bose-Einstein 
condensates, i.e., the dilute matter- wave-coherent gases formed if cooled to ultralow 
temperatures, where the critical temperatures are of order Tc ~ 100 nK • • • 1 /iK; for 
a short review of the (relatively) recent status of this rapidly developing field see 
Ref.l^. In what follows, it will be of some importance that Bose-Einstein condensates 
belong to a special class of quantum perfect fluids, so-called superfluids.!^ 

The curved space-times we have in mind are experienced by sound waves prop- 
agating on the background of a spatially and temporally inhomogeneous perfect 
fluid. Of primary importance is, flrst of all, to realize that the identification of sound 
waves propagating on an inhomogeneous background, which is itself determined by 
a solution of Euler and continuity equations, and photons propagating in a curved 
space-time, which is determined by a solution of the Einstein equations, is of a kine- 
matical nature. That is, the space-time metric is fixed externally by a background 
field obeying the laws of hydrodynamics (which is prepared by the experimental- 
ist) , and not sclf-consistently by a solution of the Einstein equations G^j/ = SttT^i, 
(where G and the speed of light are set to unity). The latter equations relate space- 
time curvature - represented by the Einstein tensor G^i, = R^^, — ^g^^R, where iJ^y 
is the Ricci tensor and R — g^'^R^n, the Ricci scalar - with the energy-momentum 
content of all other fundamental quantum fields. This energy-momentum content 
is represented by the classical quantity T^^,, which is the regularized expectation 
value of a quantum energy-momentum tensor."^ 

As a first introductory step to understand the nature of the proposed analogy, 
consider the wave equation for the velocity potential of the sound field $, which in 
a homogeneous medium at rest reads 



2 



1 9 

^TT^ + A 



$ = 0, (1) 



where Cs is the sound speed. It is a constant in space and time for such a medium. 
This equation has Lorentz invariance, that is, if we replace the speed of light by the 
speed of sound, it retains the form shown above in the new space-time co-ordinates 
obtained after Lorentz-transforming to a frame moving at a constant speed less than 
the sound speed. Just as the light field in vacuo is a proper relativistic field, sound 
is a "relativistic" field.'' The Lorentz invariance can be made more manifest by 
writing equation (^ in the form 0$ = rj^^'^di^id^^ = 0, where rj^'^ — diag(— 1, 1, 1, 1) 
is the (contravariant) flat space-time metric (we choose throughout the signature 
of the metric as specified here) , determining the fundamental light-cone structure 
of Minkowski space pi we employ the summation over equal greek indices /i, j^, • • • . 
Assuming, then, the sound speed Cg = Cs{x,t) to be local in space and time, and 
employing the curved space-time version of the 3-l-lD Laplacian n,'^one can write 
down the sound wave equation in an inhomogeneous medium in the generally co- 

''More properly, we should term this form of Lorentz invariance pseudorelativistic invariance. We 
will however use for simplicity "relativistic" as a generic term if no confusion can arise therefrom. 
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variant forni!^ 

^5^(V^g^"5,$) = 0. (2) 



Here, g = det[g^^] is the determinant of the (covariant) metric tensor. It is to 
be emphasized at this point that because the space and time derivatives 9^ are 
covariantly transforming objects in |(5J, the primary object in the condensed-matter 
identification of space-time metrics via the wave equation (jSJ is the contravariant 
metric tensor g'"'!^^In the condensed-matter understanding of analogue gravity, 
the quantities g^'^ are material- dependent coefficients. They occur in a dispersion 
relation of the form g^^kf^ki, = 0, where fc^i = i^^/cs, k) is the covariant wave vector, 
with hk the ordinary spatial momentum (or quasi- momentum in a crystal). 

The contravariant tensor components g^"^, for a perfect, irrotational liquid turn 
out to be^-^ 

' " ^ (3) 



^cc2 V -^ Ctl 



where 1 is the unit matrix and Ac a space and time dependent function, to be 
determined from the proper equations of motion for the sound field (see below). 
Inverting this expression according to g g^^ = S^ a, to obtain the covariant metric 
gfj,u, the fundamental tensor of distance reads 



v'^] —V 



where the line element is ds^ = g^^dx^dx'^ . This form of the metric has been 
derived by Unruh for an irrotational perfect fluid described by Euler and continuity 
equationsp' its properties were later on explored in more detail in particular by 
Visser.'-'J I also mention that an earlier derivation of Unruh's form of the metric 
exists, from a somewhat different perspective; it was performed by Trautman. ■ 

The conformal factor Ac in Q depends on the spatial dimension of the fluid. It 
may be unambiguously determined by considering the action of the velocity poten- 
tial fluctuations above an inhomogeneous background, identifying this action with 
the action of a minimally coupled scalar field in Z> + 1-dimensional space-time l^^^ 



S = I d^+'x- 
2g 



^$--y-V$) +ci{v^y 



= i y d^+i:Ey^g^''9^$9,$ , (5) 

where it is assumed that the compressibility 1/g of the (barotropic) fluid, g = 
d{\nnip)/dp, where p is the pressure and p the density of the fluid, is a constant. 
Using the above canonical identification, it may easily be shown that the conformal 
factor is given by Ac = (cs/g) ^ . It is mentioned here that the case of one spatial 
dimension (D = 1) is special, in that the so-called conformal invariance in two 
space-time dimensions implies that the classical equations of motion are invariant 
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(take the same form) for any space and time dependent choice of the conformal 
factor Ac- 

The Hne element ds"^ = g^^^udx^dx" gives us the distances travelled by the 
phonons in an effective space-time world in which the scalar field $ "lives" . In par- 
ticular, quasiclassical (large momentum) phonons will follow light-Vike, that is, here, 
sound-Wke geodesies in that space-time, according to ds^ — 0. Particularly notewor- 
thy is the simple fact that the constant time slices obtained by setting dt ~ in 
the line element are conformally flat, i.e. the quasiparticle world looks on constant 
time slices like the ordinary (Newtonian) lab space, with a simple Euclidean metric 
in the case of Cartesian spatial co-ordinates we display.^ All the intrinsic curvature 
of the effective space-time is therefore encoded in the metric tensor elements goo 
and goi. Together with Matt Visser, I described this curvature and its properties for 
an isotachic fluid (i.e., having a speed of sound indepen den t of space and time).'^ 
Using that phonons move on geodesies, we discovered irM^the phenomenon th at a 
vortex acts on quasiclassical phonons as an effective gravitational lens. In Ref.l^, 
using the fact that there is curvature in any spatially inhomogeneous flow (that is, 
a flow which is not a simple superposition of translational motion and rigid body 
rotation), we have shown that there exists a sonic analogue of the "warp-drive" 
in general relativity permitting superluminal, i.e., here "superphononic" motionM^ 
The point made by us is that in the acoustic curved space-times we consider, there 
is no violation of any condition on the positivity of energy necessary, which is in 
marked contrast to the original warp drives, where lo cal energy densities by ne- 
cessity must be negative to permit superluminal travel!^ This is due to the fact 
that the Einstein equations, relating curvature and energy-momentum content of 
all fields other than gravitational fields, as already mentioned in the above do not 
need to be imposed in the analogy. 

It is important to recognize that the form |(5Jl of the wave equation is valid gener- 
ally (with a possible additional scalar potential term). That is, a generally covariant 
curved space-time wave equation can be formulated not just for the velocity pertur- 
bation (sound) potential in an irrotational Euler fluid, for which we have introduced 
the effective metric concept. If the spectrum of excitation (in the local rest frame) 
is linear, uj — Cpropfc, where Cprop is the propagation speed of some collective exci- 
tation, the statement that an effective space-time metric exists is true, provided we 
only consider wave perturbations of a single scalar field which constitutes the fixed 
classical background. The argument to reach this conclusion is as follows. 

Given that the action density C is a, functional of (/> and its space-time deriva- 
tives dfj^cj), i.e. C = C[(p, 9,i</>], we expand the action to quadratic order in the fluctua- 
tions around some stationary classical background solution (Jjq of the Euler-Lagrange 
equations. For any Lagrangian of the specifled form, the wave equation for pertur- 

''This fact implies that metrics with nonvanishing spatial curvature, Uke the Kerr metric, are 
not amenable within this simple e ffect ive metric scheme, which starts from Eulcr and continuity 
equations; see for a discussion Ref.LUJ 
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bations (50 = $ above the background 0o is^ 



d, 



d^C 



^ \ 9(5^0)9(5^0) 



a.$ - 



d'^C 



-d, 



The above equation in covariant notation reads 
1 



-:d,U 



d(t)d(t) "" \ 9(9^0)90 
ion reads 

^"-9,$) - Kff(0o)$ - 0, 



* = 0. (6) 



<t>=4>o 



(7) 



where 14ff(0o) is a background-field dependent efi^ective potential (equal to the 
second term in round brackets in equation © divided by \/—g). The potential V^ff 
may, for example, in the simplest case contain an effective mass of the scalar field. 



'prop 



$. 



such that the wave equation becomes Klein-Gordon like, n$ = —rri^ 

The effective metric coefficients are, up to an (again dimension dependent) con- 
formal factor given by: 



g^ii^iM oc 



did^(b)d{d^cj,) 



(8) 



The concept of an effective space-time metric therefore applies to every system hav- 
ing a single scalar wave equation of second order in both space and time derivatives, 
corresponding to perturbations propagating on a fixed classical background, where 
this background itself determines the metric coefficients. 



2. The effective metric in Bose-Einstein condensates 

We assumed in Eq. (O that the compressibility l/g is a constant. This entails that 
the (barotropic) equation of state reads p = \gp^- We then have, in terms of the 
interaction between the particles (atoms) constituting the fluid, a contact interac- 
tion (pseudo-)potential, V{x — x') = g6{x — x'). This is indeed the case for the 
dilute atomic gases forming a Bose-Einstein condensate. Well below the transition 
temperature, they are described to good accuracy by the Gross-Pitaevskii mean 
field equation for the order parameter ^ = (^),'^ representing the expectation value 
of the quantum field operator '5: 



ih—-^{x,t) 
at 



2m 



A + Vtr^p{x,t)+g\-9{x,t)\^' 



*(a;,i). 



(9) 



The Madelung transformation reads ^ = Y^exp[i0], where p yields the condensate 
density and is the velocity potential; it allows for an interpretation of quantum 
theory in terms of hydrodynamicsM^ Namely, identifying real and imaginary parts 
on left- and right-hand sides of ©, respectively, gives us the two equations 

ft2 A^P 



^d , 


1 


^Tr(^ = 


—mv 


dV 


2 



Vtr, 



9P 



2m y/p 



P+PQ, 



(10) 



'^Observe that ('I') ^ breaks pa rticl e number conservation (the global U(l) invariance); for a 
review of the consequences see Refl^sl 
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^p + V-(H=0. (11) 

The first of tliese equations is tlie Josephson equation for tlie superfluid phase, 
which corresponds to the Bernoulli equation of classical hydrodynamics, where the 
usual velocity potential of irrotational hydrodynamics equals the superfluid phase 
(j) times h/m, such that v — K\I(j)/m. The latter equation implies that the flow 
is irrotational save on singular lines, around which the wave Junction phase <p is 
defined only modulo 27r. Therefore, circulation is quantized,'^^^ and these singular 
lines are the center lines of quantized vortices. The isothermal chemical potential 
^ (which we have chosen to incorporate the kinetic energy term ^mv'^), is aug- 
mented by the "quantum pressure term" pq = —-^{/S.^)/^, which is the one 
genuine quantum term in l|10|l . because pq ex h^ (observe that the first order in h 
may be incorporated into the velocity potential). The second equation (fTT|l is the 
continuity equation for conservation of particle number, i.e., atom number in the 
superfluid gas. The dynamics of the weakly interacting, dilute ensemble of atoms is 
thus that of a perfect Euler fluid with quantized circulation of singular vortex lines, 
which are the only vortical excitations in a superfluid. The fact that it is an Euler 
fluid is true save for regions in which the density rapidly varies and the quantum 
pressure term pq becomes relevant, which happens on scales of order the coherence 
length ^0 = ^/V^^mpo, where po is a constant asymptotic density far away from 
the density-depleted (or possibly density-enhanced) region. The quantum pressure 
becomes relevant in the depleted-density cores of quantized vortices, or at the low- 
density boundaries of the system, and is negligible outside these domains of rapidly 
varying or low density. 

The whole armoury of space-time metric description of excitations, explained 
in the last section, which is based on the Euler and continuity equations, is then 
valid for phonon-like excitations of a Bose-Einstein condensate, with the space-time 
metric (0J, as long as we are outside the core of quantized vortices,'^ and far enough 
from the boundaries of the condensate. 



3. Nonuniqueness of the quasiparticle content of a Bose-Einstein 
condensate 

Quasiparticles are the fundamental entities used to describe an interacting con- 
densed matter system in a particle picture, that is, in a suitable Fock space. 
On the microscopic level, if the elementary constituents interact by two-body 
forces, we are given a second quantized Hamiltonian operator of the form H = 
Sfc ^k<^kO,k + Vkk' cikiO'p.O'kO.k' ^ where Vkk' are the matrix elements for two-particle 
interaction in a plane wave basis, and e^ are the bare single particle energies of 
the "elementary" bosons or fermions, which are created by the operators a\, from 

"^A treatment of sound wave propagatio n in the presence of vorticity, and the corrections to Eq. 
arising therefrom, may be found in Ref.l=i| 
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the proper particle vacuum. One th en employs a unitary, i.e. operator-algebra- 
conserving Bogoliubov transformation^ to another set of quasiparticle operators 
&fc, bk' s^^ ^^- ^^'^^ below. This gives the Hamiltonian the reduced diagonal form 
^rod = J2k^(^)^k^'' + 0{b'^), where 0(6^) represents additional terms of higher 
order than quadratic, which are supposed to be small compared to the leading 
diagonalized part of the Hamiltonian, for the picture of noninteracting quasiparti- 
cles to make sense. These quasiparticles possess a (possibly spatially anisotropic) 
dispersion relation lo — uj{k) which is linear for various important classes of quasi- 
particles, Lu{k) ex k. Among these classes of collective excitations are, e.g., phonons, 
antiferromagnet i c rn agnons }^ or the excitations around the gap nodes of the p- wave 
superfluid ^e-A.E3l 

Phonons are the small momentum quanta of the sound field in solids or fluids, 
with uj{k) = Cgk for a medium at rest. Their classical equation of motion in a 
perfect fluid is the generally covariant wave equation |2Il, with the effective space- 
time metric Q. In the simplest case, for a homogeneous medium with space and 
time independent density, and a c onstant speed of sound Cg, the quantized velocity 
potential of the phonons reada^ 

~iCgkt+ik-x L^ 1^ iCskt — ik-x 




where V is the quantization volume of the system, and po the constant background 
density. That is, the quantum excitation field in a homogeneous medium may gener- 
ally be decomposed into plane waves, with the appropriate frequencies as a function 
of momentum stemming from the dispersion relation, here lo — Cgk. Note, in partic- 
ular, that for this spatially and temporally homogeneous case, the statement that 
we observe positive frequency (energy) with respect to the laboratory time interval 
dt is unique, that is, it can be made independent of time and space. In an inhomo- 
geneous fluid, this is (generally) no longer the case, and the notion of an excitation 
having positive energy may depend on where the detector is located in the fluid, if 
it is at rest relative to the fluid or moves, and what its natural time interval is. The 
latter may be different from that of the laboratory, due to the particular way the 
detector couples to the fluid, see sectional 

3.1. Operator basis dependence of quasiparticle content 

The number of particles assigned to the quantum field $ is unique with respect 
to a certain given state |6) of the quantum field, nk{&) = (6|aj,afe|9), provided 
we decompose <& into modes associated to the operators a^ and their Hermitian 
conjugates. The number of particles is, in particular, zero for the vacuum state 
with respect to the field operator dk, defined by dk\0) — 0. However, it need not be 
zero with respect to another set of quasiparticle operators bk , which has a different 
vacuum |0). 

To demonstrate the basis dependence of quasiparticle number, we use a general 
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Bogoliubov transformation for bosons, the class of (quasi-)particles we consider, of 
the general form 



a 

k 



y^ aikbk + Pikb\ ' (1"^) 



where k represents a set of quantum numbers, not necessarily the plane waves used 
in Eq. H12|l . The coefficients in the Bogoliubov transformation must fulfill certain 
conditions for the transformation to be unitary, i.e., to preserve the bosonic com- 
mutation relations for the new operators, [6i,6^] = 5ik, [6i,6fc] = 0, [^1,^1] = 0. 
As a consequence of this defining unitary character, the following conditions on the 
transformation coefficients must hold: 

^ aika*jf, - (3ikf3jk = Sij, ^ aikPjk - Pikajk = 0. (14) 

k k 

By using the transformation H13() . it is straightforwardly shown that the number 

of dfc particles in |0) is given by (0|a^afc|0) = J2k' l/^fc'fcP- The old operator Uk 

does not annihilate the new vacuum |0) (and vice versa), and what looks empty in 

one quasiparticle vacuum may be full of quasiparticles in another. Related to this 

(formal) operator-basis dependence is the fact that the actually detected number of 

particles is strongly observer dependent, as opposed to the formally defined quantity 

'^fc(0), which refers to one particular quasiparticle state |6). The detected number 

of particles depends, in particular, on how the detector actually couples to the 

field <I> whose quanta it measures. Various couplings of the detector, for example to 

different powers of the fluid density, will influence the quasiparticle basis in which 

the detector measures, and thus the quasiparticle number detected. 

Now, the salient point is that because the phonon field $ is a relativistic quantum 

field, we will be able to map the observer dependence just described, which is general 

and holds for any sort of quasiparticles, to the observer depend ence ex perienced by 

proper relativistic quantum fields in curved or flat space-time! ' ' The observer 

dependence is of kincmatical origin, i.e., it originates in the fact that the relativistic 

quantum field propagates in a space-time of Lorcntzian signature, cf. the discussion 

on 
of the Hawking radiation analogue phenomenon in Ref. , and is therefore fully 

within the capabilities of our proposed analogy. 

3.2. Scaling ansatz in expanding Bose-Einstein condensates 

To model the quasiparticle analogue of expanding universes, we will make use of 
expanding Bose-Einstein condensates, which are produced by reducing the trapping 
potential strength, i.e. the harmonic trapping frequency with time, or by increasing 
the interaction coupling constant. Firstly, it is thus appropriate to describe the evo- 
lution of density and velocity distribution in t he exp anding gas by discussing the 
so-called scaling procedure established in Refs.l I I The scaling procedure intro- 
duces a set of generally three scaling variables, bi = bi{t), which are a function of 
time only. Using these scaling variables, one writes for the (Cartesian) co-ordinate 
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vector components xi,i — Xi/bi] for the scaled co-ordinate vector we use the short- 
hand notation Xb = 'Y^^eiXi/bi. It may then be shown that the evolution of the 
Bose-Einstein-condensed gas cloud is described, starting from the initial density 
and velocity potential distributions p = pinit{x,t = 0), ^ = 4>iait{x,t = 0), by the 
following density and velocity distributions (from here on we take h — ra — i)lsisilMI 

P(^,t) ^ ^i^, (15) 

v^V(t){x,t) => V = V(j) = ^-^Xi+y(i>{xb,t). (16) 

This is true provided the scahng parameters 6^ fulfill the equationg^^l^ 

&,+c.f(i)6. = 4^^, (17) 

.9(0) Vhi 

where wqi are the initial trapping frequencies and the dimensionless "scaling vol- 
ume" reads V — J^,- bi . The dots are time derivatives with respect to laboratory 
time t here and in what follows.'^ We have taken into account in the above equation 
that the particle interact ion can be varied in lab time, g = g{t), by means of a suit- 
able Feshbach resonance) I ' g{0) is the initial coupling constant. We designate 
"scaling basis" quantities with a tilde. For example, the (stationary) initial density 
distribution ^11111(3^) gives us Pinitixb) if we replace x —> x^. The scaling evolution 
is exact in a Thomas-Fermi approximation which neglects the quantum pressure 
term pq and the kinetic energy ^mv^ . Since Bose-Einstein condensates were exper- 
imentally created, the scaling solution has routinely been employed to interpret the 
time-of-flight pictures with which they are visualized as well as analyzed.l^l 
We define a "scaling time" variable by 

dTs 5(t)/.9(0) 



dt V 

and the Ts dependent scaling functions Fi by 



(18) 



f'i 9{rs) bf dTs 
In terms of these quantities, the effective second order action for the scaling basis 
fluctuations of the phase of the superfluid order parameter 5(l){xb,Ts) = ^{xb,Ts) 
takes on the particularly simple diagonal forni^ 



S'(2) = f dTsd^Xb 



25(0) 



—(5^) +d^F,(yb 



(20) 



where the scaling speed of sound c — ^yg{Q)pinit{xb)■ Because of the fact that this 
action does not mix spatial and temporal derivatives, the resulting line element 

"Note that by special convention we generally do not sum over latin indices i,j, ■ ■ ■, but only over 
greek indices fj,,u,- ■ ■. 
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10 

in the scaling variables, according to the identification displayed in ((SJ is diagonal 
(does not possess goi terms), and reads 



ds^ = g^^udx^^dx" = —-^F^FyF, [-i?dT^, + FrHxl] . (21) 

We now consider for simplicity the isotropic case, hi = h, which implies 

i^, = ^ = 6^-23(0) _ ^22) 

The generalization of the mode expansion in Eq. p2|) to inhomogeneous expanding 
Bose-Einstein condensates then takes, in this isotropic case, the forni^ 



5(0) 



$(a=,,T,) = ^Jii^0„(a;b) kxn{Ts)+hixl{rs) . (23) 

n V ^^^" 

The functions 0„(a;b) are the stationary solutions of the Gross-Pitaevskii equations 
for excitations above the initial ground state, designated by the (set of) quantum 
numbers n with initial energies e„. The initial Thomas-Fermi quantization volume 
is V; in the hard-walled cubic box limit the modes are plane waves, 0„(a;b) -^ 
exp[ifc • Xb]. 

The tem por al mode functions Xn{Ts) satisfy the second order ordinary differen- 
tial equation^ 

^Xn + F{Ts)elxn = 0. (24) 

The case when _F is a constant (unity) is particular. In this case, the quantum state of 
the quasiparticle excitations remains unchanged for increasing Ts, and a given initial 
quasiparticle vacuum, in the scaling basis with the associated quasiparticle operators 
defined by (|23|l . remains empty forever. That is, no quasiparticles are created in that 
basis, although the superfluid may be in a highly nonstationary motional state, 
obtained by changing the trapping lo — Lu(t) rapidly with time. However, a detector 
which measures in a quasiparticle basis different from the scaling basis, for example 
due to its particular coupling to the superfluid, may still detect that quasiparticles 
are "created" . We will come back to this possibility in section 01 below, when we 
discuss the purely choice-of-observer related phenomenon of a thermal state in a 
quasiparticle basis belonging to one particular space-time, the de Sitter space-time 
of cosmology. 



3.3. "Cosmological" quasiparticle production 

Consider now the general case that the scaling function F{ts) is a function of scaling 
time Tj. The fact that F depends on time implies that the statement "the excita- 
tion is of positive frequency" (a particle) or of "negative frequency" (antiparticle) 
for a given propagating wave cannot be held up for all times Tj. This frequency 
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11 



mixing implies that quasiparticles are created from the quasiparticle vacuum, be- 
cause an initiaUy empty scaHng basis vacuum state does not remain empty during 
the evolution of the system, i.e. initially 6„|0(rs = 0)) = 0, but at a later stage 

KMts)) ^ 0. 

The fact that annihilation and creation operator parts of the initial vacuum are 
mixed, as a consequence of H13|l . is physically due to the fact that quasiparticles 
are scattered within the course of time at (time dependent) effective potentials. 
Physically, the Tg dependence of F furnishes such an effective potential, at which 
excitations are scattered from negative to positive frequency and vice versa: The 
equation 1)24(1 is formally equivalent to scattering of a non-relativistic particle with 
energy e„ by a potential in r^. space. 



V{rs)^el{l-F{rs)). 
At large Tg, the WKB scattering solution of 124|l therefore readsM^ 
1 



Xri 



Fi/4 



a„exp 



len / driy^Fi^ 



+ Pn exp 



ze„ / dri^/Wi^) 



(25) 



(26) 



where the scattering amplitudes are related via the particle flux conservation con- 
dition |a„p — |/3„p = 1. The quantity Nn = |/3„p can be interpreted as the number 
of scaling basis quasiparticles created from the initially empty scaling vacuum, due 
to the time dependent scattering of excitations moving in the nonstationary con- 
densate. 

In the WKB approximation, the amplitudes a„ and /3„ are connected in a simple 
way: 



/3„ = exp 






where the inverse temperature I/Tq is given by the integral 

1 



^ 



Fdr, 



(27) 



(28) 



and C is the contour in the complex r^ -plane enclosing the closest to the real axis 
singular point of the function a/F(ts)!^ This gives us the number of particles 
created in the mode n, using the flux conservation condition |a„p — |/3„p = 1, 



A^„ 



\Pn 



(29) 



exp[e„/ro] - 1' 

The distribution of the created quasiparticles follows a thermal bosonic distribution 
(Planck spectrum), at a temperature Tg. The adiabatic evolution of trapped gases 
hence leads to "cosmological" quasiparticle creation with thermal occupation num- 
bers in the scaling basis. The temperature Tq occurring in the Planck distribution 
above depends on the details of the scaling evolution, i.e., on the specific superfluid 
dynamics imposed by the solution of Eqs. H17|) . that is, To is a functional of the 
temporal evolution Wj = uji(t) . "^^^ 
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We use the term "cosmological" in context with the plain condensed-matter fact 
that quasiparticles are created in the scahng basis. We now justify this by comparing 
our effective Bose-Einstein condensate metric H21|l to hne elements which constitute 
cosmological solutions of the Einstein equations. For example, let dr^ = dt by 
properly adjusting g = g(t), thus choosing the coupling constant's time evolution to 
be given by g{t) = g{0)V{t), cf. Eq. l(TF|l . We then obtain that l(?T)) equals (up to the 
conformal factor c/{g(0)V)) an anisotropic version of the spatially flat Friedmann- 
Robertson- Walker (FRW) universe?^ 



ds' 



5(0)V 



-c^dt" +Y.^1d4i 



(30) 



In the standard spatially isotropic form of the FRW metric, all bi = b are equal. 

To obtain the exact equivalence to a spatially flat FRW metric, we have to 
assume in addition that c is spatially independent, which is fulfilled close to the 
center of the gas cloud, where the parabolic density profile is approximately flat, 
and c is essentially a constant. In the spatially isotropic case, 

H^l (31) 

is the Hubble "constant" , which obviously is a constant only if the gas expands ex- 
ponentially in laboratory time, b oc exp[Ht\, just as we need exponentially rap id ex:- 
pansion for a constant Hubble parameter in inflationary cosmological models! I '^ ^ 
We therefore come to the remarkable conclusion that the co-ordinate scaling 
factor b of the Bose-Einstein condensate quasiparticle universe, occurring in the 
equations of motion of a nonrelativistic condensed matter system, has exactly the 
same meaning, in the quasiparticle world, as the cosmological scale factor of the 
Universe proper. 



4. Gibbons-Hawking effect in de Sitter space-time 

We now treat the case that the scaling factors Fi in l|19(l are constants, i.e., do not 
depend on scaling time. Therefore, following Eqs. (|23|l and H24|l . no scaling basis 
quasiparticles are created through negative and positive frequency mixing. The 
superfluid can still be in highly nonstationary motion, though: The time evolution 
is according to iflt^ prescribed by 

9W Ih^iOk g(t) 

cf. Eq. H19|l . where the Ci are constants. However, although the superfluid is in 
motion, no dissipation through intrinsic quasiparticle creation takes place, because 
there exists the Fock space "scaling" basis, in which no quasiparticles are created 
from the scaling vacuum, and the energy of that particular superfluid vacuum is 



February 2, 2008 8:27 WSPC/INSTRUCTION FILE FINALMPLAreview 



13 

conserved. A particular instance is the isotropic 2D case, bi — b2 — b, where for 
constant g{t) — g{0) the condition on F being constant is fulfilled/ 

In de Sitter space-time space is empty and flat (that is, the constant time slices 
are Euclidean space), and all of the curvature of space-time is encoded in a nonvan- 
ishing cosmological constant A. For obvious reasons, the de Sitter space-ti me is v ery 
popular in the quantum field theoretical treatment of cosmological theories) II be- 
cause it highlights the crucial cosmological role the energy density of all conceivable 
quantum fields taken together might play: The vacuum energy density may con- 
stitute the dominant effective source term for space-time curvature in the Einstein 
equations. 

The Gibbons-Hawking effect for geodesic observer^^ in such a de Si tter space- 
time is the curved space-time analogue of the Unruh-Davies effect! 1 ^ 1 The latter 
consists in the fact that a constantly accelerated detector moving in the fiat (pur- 
portedly "empty") space-time vacuum, responds as if it were placed in a ther- 
mal bath of (quasi-)particles with temperature proportional to its acceleration. 
Observer-related phenomena are at the heart of quantum field theory on nontrivial, 
and generally curved, space-time backgrounds. They tell us, in particular, that the 
particle content of a given quantum field depends on the (motional) state of the 
detection apparatus, which is verifying that there are particles by its "clicks". More 
technically speaking, the dependence of the particle content of quantum fields in 
curved space-time is rooted in the non-uniqueness of canonical field quantization 
in Riemannian spaces^ It is of fundamental importance to make these observer- 
dependent effects measurable, because such a measurement constitutes, inter alia, 
a consistency check for an all-important concept of standard quantum field theory, 
namely that the quantization of a given field is carried out on a fixed space-time 
background. 

That the experimental verification of the observer dependence is exceedingly 
difficult with light becomes readily apparent if we calculate the Unruh-Davies tem- 
perature: The result is that it equals Tunruh — [ft/(27r/cBCL)]a — iK x a[10^°g©] , 
where a is the acceleration of the detector in Minkowski space (g® is the gravity 
acceleration on the surface of the Earth), and cl the speed of light. The huge ac- 
celerations needed to obtain measurable values of Tunruh make it obvious that an 
observation of the effect with light (photons) is decidedly a less than trivial un- 
dertaking. Although proposals for a measurement with ultrai ntense short pulses of 
electromagnetic radiation have been put forward in, e.g., Refs! '"^ 1 '^ I, it is less than 
obvious how the thermal spectra associated to the effect, which still furnish tiny 
contributions to the total energy, should be discernible from the background domi- 
nated by the ultraintense lasers used to create large accelerations of the elementary 
particles contained in the plasma. 

*Cf. Ref.lss^ where the fact of superfiuid vacu um energy conservation is explained from a differ- 
ent (SO(2,l) symmetry) perspective, and Ref .ti2l for a quahty factor measurement of breathing 
(monopolc) oscillations in a cylindrical geometry. 
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4.1. An isotropic de Sitter universe in a harmonic trap 

We first discuss the simplest case of an isotropically expanding gas in a liarmonic 
trap. We will see that this case also serves a pedagogical purpose, because it forces us 
to distinguish between purely observer-related phenomena with thermal spectrum, 
which by definition all have Fi^Tg) = 1, and "cosmological" particle creation with a 
thermal spectrum, as defined in section rTlH for which the scaling functions Fi{Ts) 
in H19|) depend explicitly on scaling time. 

One can create de Sitter universes in an expanding gas by letting it expand expo- 
nentially, h{t) ex exp[iJi] where H is the Hubble constant of cosmological expansion, 
equalling H = A/D for the de Sitter universe discussed, where D is the spatial di- 
mension. The de Sitter metric i n it s standard form, using as its time co-ordinate 
the "cosmological" time Tc, reada^ 



where cq = c(0) is the central scaling (i.e., initial) speed of sound. In the present case 
of exponential expansion, the "cosmological" time interval occurring in the metric 
above is equal to both the laboratory and the scaling time interval, dTc = drg = dt; 
cf. Eq. ^ with bi = b = exp[iJrc]. 

However, using the isotropic harmonic expansion setup, one has to face the dif- 
ficulty that in order to give the FRW metric (|30|l near the center of the trap (i.e., 
close to Xb — 0) the de Sitter form, one has to increase exponentially the interaction 
with laboratory time: Because of (|31|l . g{t) (x b^ (x exp[DHt] = exp[coAi]. Though 
the central density also decreases exponentially (like b~^^^), the exponential in- 
crease of the coupling "constant" incurs strong three-body recombination losses,!^ 
whose total rate (in the dilute gas case and in three spatial dimensions) is propor- 
tional to g^p^ oc 6^ = exp[9coAt]. Therefore, within a short time of order 1/H, the 
Bose-Einstein-condensed gas of interacting single atoms will simply no longer be in 
existence, because the atoms rapidly form bound states. Such an experiment will 
leave no time to measure phenomena which depend on the fact that an equilibrium 
is established; in particular, the thermal equilibrium for the occupation numbers in 
the de Sitter quasiparticle basis will not be established on such a short time scale. 

Even more importantly, though one obtains indeed a thermal spectrum in the 
de Sitter quasiparticle basis corresponding to the metric (|33|l . F{ts) = F{t) = 
1/6^ = exp[— 2_ffi] depends on time, i.e., it is not a constant in the quasiparticle 
basis corresponding to the mode functions Xn oc exp[— ie„rs]. The thermal spectrum 
obtained therefore is, by our physical definition, the thermal "cosmological" quasi- 
particle creation discussed in the previous section. It is not what we want to observe 
in our effective de Sitter space-time, namely the purely choice-of-observer related 
phenomenon Gibbons-Hawking effect, for which no actual quasiparticle "creation" 
in the scaling basis should take place. 
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4.2. The 1 + lD de Sitter universe in a cigar-shaped cloud 

To circumvent the problem that the interaction couphng needs to be increased 
exponentially with time to obtain a de Sitter universe in 2+1 or 3+1 isotropic 
space-time dimensions in an isotropic harmonic trap, I and Petr Fedichev have 
developed the model of a 1+lD de Sitter universe. This 1+lD toy mode l can be 
realized in a strongly anisotropic, cigar-shaped Bose-Einstein condensate) I I cf. 
Fig-D III particular, in the proposed experimental setup, no time variation of the 
coupling constant at all is necessary, which is thus a true "constant" also in time. 
The analysis of the excitation modes in a strongly anisotro pic, elongated Bose- 
Einstein condensate is based on the adiabatic separation ansatz^ 

$(r,z,t)=^0„(r)x„(z,t), (34) 

n 

where (/)n(r) is the radial wavefunction characterized by the quantum number n 
(only zero angular momentum modes are considered here) . The above ansatz incor- 
porates the fact that for strongly elongated traps, i.e., traps for which uj^ <C ujj_, the 
dynamics of the condensate motion separates into a fast radial motion and a slow 
axial motion, which are essentially independent. The x„(z,t) are the mode func- 
tions for travelling wave solutions in the z direction (plane waves for a condensate 
at rest read Xn oc exp[— je„.fci + fcz]). The radial motion is assumed to be "stiff" such 
that the radial part is effectively time independent, because the radial time scale 
for adjustment of the density distribution after a perturbation is much less than 
the axial oscillation time scales of interest. The ansatz H34(l works independent from 
the ratio of healing length and radial size of the superfluid cigar. In the limit that 
the healing length is much less than the radial size, Thomas-Fermi wave functions 
are used, in the opposite limit, a Gaussian ansatz for the radial part of the wave 
function (pni^) is appropriate. 

The squared oscillation spectrum of the ci gar- shaped condensate cloud reads 
^n k — ^0^^ + 2uj\^n{n + 1), where Cq = ^/J^/2<^ A de Sitter space-time is then 



b± (X \/t 




Atomic 
Quantum Dot 



b (X t 



Fig. 1. Expansion of a cigar-shaped Bose-Einstein condensate, with scale factor b Hnear in lab 
time along the axial, and with •/! in the radial direction, scaled by b±. The stationary horizon 
surfaces are located at ±zh, respectively. The thick dark lines represent lasers creating an optical 
potential well in the center of the harmonic trap, which hosts the Atomic Quantum Dot, cf. Fig.l2l 
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obtained from the effective action for the phase fluctuations of the phonon {n = 0) 
modes near the center of the trap. The corresponding I+ID action is o btain ed after 
integrating out the transverse, strongly confined directions, and reads ! I I 



So ~ dtdz 



^j_jij_ 



liX.~v.d) +^(5.Xof 



25 
^ i I d^+ia;V^g^''a,,xoa.xo . (35) 

The scaling parameters h in the axial (wz) and h±^ in the perpendicular {ijJ±_) direc- 
tions are functions of time, such that the action fulfills the identification with the 
action of a scalar field minimally coupled to gravity, as indicated in the second line 
above, where the metric coefficients are the g^^ in Eq. I^J. 

The identification of the above phase-fluctuations action with the action of a 
scalar minimally coupled to gravity works, that is, the two actions in the first and 
second lines of (|35|l are consistent, if we impose the consistency condition that 

9 Co 

b I „ [tt 1 



where Z is a renormalization factor according to xo = Zxoj with xo the renormalized 
wave function, and pm the initial central density. The factor Z does not influence the 
(classical) equation of motion SS/Sxo = (it simply drops out) , but does influence 
the response of a detector. We will come back to this point in section 14.31 below, 
when we discuss the explicit dependence on Z of the equilibration time scale for the 
stationary detec tor state considered, see Eq. (^HJ. 

In Refs! I I, we were using an alternative form of the 1+lD de Sitter metric J33II. 
This alternative form is the one used in the original Gibbons-Hawking paper.l^ It 
reads (we leave out the conformal factor) 

ds^ = -cl (1 - Az2) dT^ + (1 - Az2) -' dz^ . (37) 

The time interval dr in the above metric is not the "cosmological" time interval dTc 
in the version of the metric displayed in (|33|l . The two metrics may be transformed 
into each other using the following co-ordinate transformations: 

exp[— 2A/Acor] — exp[— 2VAcoTc] — Az^, 

z = Zbexp[AcorJ, (38) 

giving us T = t(tc, Zf,), z = z{tc, Zf,), which transforms (|37|l into (|33|l . The advantage 
of the form H37|l is that it is plain in this form that the de Sitter space-time has an 
event horizon, located in our 1+lD case at the constant values z = zh = ±A~^^^. 
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The quantity Az^ in the de Sitter metric H37|) must be independent of time. This 
leads us to the requirement 

Va.= ^ = -^z=^z (39) 

c(0 Vbco Co 

relating the dynamical parameters of our problem to each other (c(i) is the instan- 
taneous sound velocity at the center of the cloud). These relations imply that b ~ 
const., and thus that b oc t. Then, the cosmological constant A becomes independent 
of time, as is necessary for an analogue de Sitter space-time to be established. 

The experimental procedure now is determined to be as follows: Prepare 
a Thomas-Fermi (i.e., sufhciently large), strongly anisotropically trapped, cigar- 
shaped condensate. Let it expand in the axial direction linear in lab time by changing 
the trapping according to the scaling equations (|17|1 . such that 6 ex t, and simul- 
taneously expand in the perpendicular direction with the square root of lab time, 
6_L oc V<, such that B = b^/vb in (|36|l is a constant. Then, a detector "tuned" to the 
de Sitter space-time (jSZl, i.e., which works in the de Sitter quasiparticle basis, will 
measure a thermal quasiparticle spectrum, with the de Sitter temperature!^"^^-^ 

rds-|^VA=f6. (40) 

The fact that a thermal spectrum is obtained can be directly derived from the 
equations of superfluid hydrodynamics, as expounded in Refs. ' ' '^^'- - ^ ! That is, it is 
not simply postulated due to the (kinematical) analogy with quantum field theory 
in de Sitter space-time, but is a result of the quantized hydrodynamic equations 
determining the evolution of the quasiparticle content of the sound field in the de 
Sitter basis. 

The relation between de Sitter time r and the laboratory time is fixed by dr = 
dt/{^/bb±^) — dt/{b{t)B) = dt/{btB). [Note that dr and the scaling time interval dr^ 
defined in Eq. (|18|l differ; drg — dt/B^b^ = dr/Bb.] The transformation law between 
t and the de Sitter time r (on a constant z detector, such that dt = dt), is given by 

— =exp[S6r], (41) 

to 

where the unit of lab time to ^ uj7 is set by the initial conditions for the scaling 

variables b and b±. 

It is important to recognize that an effective exponential "acceleration" of the 

oscillation frequencies, either because of the exponential dependence of laboratory 

time on scaling time, represented by 1)41(1 . or coming from the WKB approximation 

for the scattering amplitudes, and the exponentially small mixing of positive and 

negative frequency parts resulting therefrom, Eq. I|27|) . is sufficient for the thermal 

spectrum to be obtained. In particular, though we have discussed a space-time which 

possesses a horizon, no pre-existing horizons in the parent space-time are necessary 

per se for thermal occu pation number distributions to be established, as has been 

pointed out in Ref.l^. An explicit example is the Unruh-Davies effect, where this 

parent, global space-time is simply Minkowski space. 
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4.3. Detecting the thermal de Sitter spectrum 

To detect the Gibbons-Hawking effect in de Sitter space-time, one has to set up a 
detector which measures frequencies in units of the inverse de Sitter time r, rather 
than in units of the inverse laboratory time t; this corresponds to measuring in the 
proper de Sitter quasiparticle vacuum, where, in particular, positive and negative 
frequency are defined with respect to t. Only then does one detect quasiparticles 
which are defined with respect to the de Sitter quasiparticle basis, and refer to a 
vacuum corresponding to exactly that space-time. 

In Refs. , I and Petr Fedichev have provided such a detector. We have shown 
that a "de Sitter basis" detector is realized by an "Atomic Quantum Dot" (AQD) 
(for a detailed exposition of AQD properties cf. Ref.l^^. The AQD can be imple- 
mented in a Bose gas of atoms possessing two hyperfine ground states a and f3; the 
level scheme is represented in Fig.|21 The atoms in the state a represent the ex- 
panding Bose-Einstein condensate, and are used to model the expanding de Sitter 
universe. The AQD itself is formed by trapping atoms in the state /3 in a tightly 
confining optical potential Vopt created by a laser at the center of the cloud. The 
interaction of atoms in the two internal levels is described by a set of coupling pa- 
rameters gcd = 47racd (c, d = {a,P}), where ad are the s-wave scattering lengths 
characterizing short-range intra- and inter-species collisions; gaa = 9, o,aa = Q-s, and 
gai3 = g. The on-site repulsion between the atoms (3 in the dot is given by the en- 
ergy level spacing U ^ gpp/l^ between states with a occupation difference of one /3 
atom, where I is the characteristic size of the ground state wavefunction of atoms (3 
localized in Vopt. We consider the so-called coUisional blockade limit of large C/ > 0, 
where only one atom of type j3 can be trapped in the dot. This limit assumes that U 
is much larger than all other relevant frequency scales in the dynamics of both the 
AQD and the expanding superfluid, and corresponds to a large "Coulomb blockade 
gap" in electronic quantum dots.'""' As a result of these assumptions, the collective 
co-ordinate of the AQD is modeled by a pseudo-spin-1/2 degree of freedom fy, with 
spin-up/spin-down state corresponding to occupation of the AQD by a single atom 
or no atom in the hyperfine state /3. A Rabi laser of frequency $7, with a detuning 
A from resonance between the two hyperfine levels a and /3, couples atoms of the 
hyperfine species a, constituting the expanding cigar-shaped superfluid, into the 
AQD. 

The detector Lagrangian reads (Ref.l^, cf. the Hamiltonian formulation in 



Ref.El): 



^AQD = i ( j^n*\ V - n^po{o,t)p{f] + 77*) 

-A + {g - g)poiO,t) + g6p + -c 



fj*fj. (42) 



The fact that the detector has the de Sitter basis as its "natural" quasiparticle basis, 
and therefore measures in de Sitter time, is due to the fact that the term linear 



February 2, 2008 8:27 WSPC/INSTRUCTION FILE FINALMPLAreview 



19 




Fig. 2. Level scheme of the "Atomic Quantum Dot", which is embedded in the supcrfluid cigar, 
and created by an optical well for atoms of a hyperfine species different from that of the condensate. 
Double occupation of the dot is prevented by a coUisional blockade mechanism. 



in the detector co-ordinate fj in the Lagrangian couples in a certain manner to the 
superfluid cigar in which it is embedded. Namely, the laser with frequency fl, causing 
transitions between the two hyperfine levels a and f3, couples to the square root of 
the central mean-field particle density. This particular coupling (represented by the 
second term in the first line of the above Lagrangian) is what we need, because 
VbB 



\/po(0,t) — ,/Jhn/hB = y^p^/{bBt) = y^Jh^dr/dt. The fact that the coupling 
coefficient is proportional to dr/dt is required to establish that the detector can 
work as a de Sitter detector, because it transforms the detector equations, which 
are obviously a priori in laboratory time, into equations in de Sitter time; see the 
equations (43) for the temporal evolution of the detector level occupations below. 
Adjusting the detuning A properly, such that A(t) = {g — g)po{0,t) — {g — 
g)pm/ib^B^t^), the first and the second term in the square brackets in the second 
line of (|42() cancel. One then obtains a simple set of coupled equations for the 
occupation amplitudes of the state tp = ipplf]) + ilja\c() of the AQD: 






. dlpa _ Wo 
' dT ~ 2 



i^i3, 



(43) 



where dr is the de Sitter time interval. Were it not for the density oscillations in the 
cigar, represented by the potential 5V , the above equations 1)43(1 would represent a 
simple two- level system, with frequency splitting wq = '^^yPrrJ^- The density os- 
cillations contained in the perturbation operator SV{t) = [g — g) Bb{T)Sp{T) cause 
transitions between the two undisturbed detector eigenstates |±) = (|a) ± |/3))/v2 
of the two-level system, which are separated by the energy uq. The density pertur- 
bations in the expanding host superfluid lead to a damping of the Rabi oscillations 
with frequency luq between these two states. This constitutes the effect of the de 
Sitter thermal bath to be observed, where the damping happens on the time scale 
displayed in H48|l below. 
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The response of the detector, that is, the transition rates between the detector 
states, can be calculated by evaluating a response function^ which makes use of 
the expectation value of the product of two 6(r)/5(r) operators. The probability per 
unit time for excitation (P-f, transition from |+) to |— ) ) and de-excitation (P_, 
transition from |— ) to |+)) of the detector takes the form! ^ 1 ^0 1 

^ = lim 1 / / drdr' {6V{T)6V{T'))e^'^'>^-'-^'^ 
dT r^tx, T J J 

= lim ^ ^^^^^ f f drdr' {b{T)5p(T)b{T')5p{T'))e^"^°'^^-''\ (44) 
r^oo T J J 

The second-quantized solution of the hydrodynamic equations for the density fluc- 
tuations above the superfluid ground state in the expanding cigar-shaped Bose- 
Einstein condensate read J ' ' 



'^^-Z.^i^^i^p^al'^'^-^'P 






H.c. (45) 



Using this solution, and inserting into (|44|l , we have show n^ l I that at late times 
T the transition probabilities per unit de Sitter, i.e., per unit detector time satisfy 
detailed balance conditions. They correspond to thermodynamic equilibrium, at the 
temperature Tds displayed in (|in|) : 

dP+/dT ^ ne , . . 

dP-/dT l + ne' ^ ' 

where the Bose-Planck distribution function takes the form familiar from thcrmo- 
statistics: 

exp[LJo/rdsJ - IJ 
The frequency wq oc fi can be varied by changing the undressed Rabi frequency 
n, varying the intensity of the Rabi laser. The detector thus has a changeable 
and therefore tunable frequency standard, which can be adjusted to scan the above 
distribution function for a given de Sitter temperature T^q. 

From the relation 1461) . we come to the remarkable conclusion that a properly 
designed detector can "see" a thermal equilibrium distribution in its quasiparticle 
basis, though it is embedded in a highly nonstationary system with respect to the 
laboratory frame. The rapidly expanding Bose-Einstein condensate represents this 
highly nonstationary system, which hosts the de Sitter quasiparticle detector AQD. I 
stress here again that H46(l is an exact result obtained by quantizing hydrodynamic 
fluctuations in a nonstationary superfluid, and not just concluded from a mere 
comparison of the phonon dynamics in our expanding superfluid with the quantum 
field theory of photons in de Sitter space-time. 

The equilibration time scale of the detector, and thus the time scale on which the 

Rabi oscillations between the detector states are damped out, is set by the detector 

frequency standard (the level spacing) wq, and by the renormalization factor ZM^ 

Tcquii = Z-^uj^^ (X (p™a3)-i/2 {^/oj^f uj^\ (48) 
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The renornialization factor Z contained in (|36|l determines the equihbration rapidity 
because it physically expresses the strength of detector-field coupling. It is related 
to the initial diluteness parameter -Dp(O) = (pma^)^^^ of the Bose-Einstein conden- 
sate and to the ratio ^/w^, which determines inasmuch the system is effectively 
one-dimensional, Z^ ex Dp{0){iL)^/ jif' . To obtain sufRciently fast equilibration, the 
condensate thus has to be initially not too dilute as well as close to the quasi-lD 
regime, for which the transverse harmonic oscillator energy scale is of order the en- 
ergy per particle, /i ^ uj±^. These two conditions have another important implication. 
The ratio of the instantaneous coherence length ^c(i) = (87rpo(0, t)as)^^/^ oc t and 
the location of the horizons z — z-^ — ±A^^", which are stationary in the present 
setup, has to remain less than unity within the equilibration time scale. ^ If this is 
not the case, the coherence length, which plays the role of the analogue Planck scale 
(which is time dependent here), exceeds the length scale of the horizon at equilibra- 
tion, and the concept of "relativistic" phonons propagating on a fixed curved space- 
time background with local Lorentz symmetry becomes invalid. The ratio S,c(t)/ zh 
at the lab equilibration time scale t = tcquii = toex.p[2Tr{Tds/^o) (m/^-l)" -Dp ^(0)] 
following from the de Sitter equilibration time in Eq. I|48f) . expressed in parameters 
relevant to the experiment, is given by 



Cc(^cquil) ""^O^dS 

— ^— i — - — ^i^ exp 



ZTT- 



^0 V^-L/ -Dp(O) 



(49) 



We see that this ratio changes exponentially with both the initial diluteness pa- 
rameter £'p(0) and the quasi- ID parameter ij,/luj_. In most currently realized Bose- 
Einstein condensates,!^ the diluteness parameter Dp ^ 10^^. Here, we initially need 
Dp(0) ~ 0(1) to have the condition ^c(tcquii)/^H < 1 fulfilled, assuming a reason- 
ably large value of the de Sitter temperature T^g. Though the condensate has to 
be initially quite dense, it is to be stressed that the central density decays like t~^ 
during expansion. Therefore, the rate of three-body recombination losses quickly 
decreases during the expansion of the gas, and the initially relatively dense Bose- 
Einstein condensate, which would rapidly decay if left with a Dp close to unity, can 
live sufficiently long, the total rate of three-body losses decreasing like Po{0, t) ex t^^. 

5. Summary 

The primary statement to be drawn from the present article is that phonons, i.e., 
low-energy linear-dispersion quasiparticles, moving in a spatially and temporally 
inhomogeneous Bose-Einstein-condensed superfluid gas, are equivalent to photons, 
the quanta of the electromagnetic field, moving on geodesies in a given curved 
space-time. We have explored the classical as well as the quantum aspects of this 
statement. 

SI thank C. Zimmermann for a pertinent question during a talk given by me at Tiibingen, leading 
to this observation. 
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On the classical side, the analogy helps to provide us with a simple general 
means to study quasiparticle propagation in an inhomogeneous medium in motion. 
An example of s uch an application is the gravitational lensing effect exerted by a 
superfluid vortex.l^ On the quantum field theoretical side, we can access within the 
analogy phenomena which are extremely difficult if not impossible to access with 
light. One of the fundamentals of quantum field theory, the fact that the particle 
content of a quantum field depends on the observer, can thus be experimentally 
verified for the first time. The basic reason that the phenomena in question are 
(comparatively) easy to simulate in a condensed matter system is that the energy 
and temperature scales, under which they occur, relative to the typical energy scales 
of the system, can be changed at will by the experimentalist in a very controlled 
manner. More particularly, the temperature of the thermal spectrum of phonons to 
be measured in the Gibbons-Hawking effect can be made relatively large compared 
to the axial phonon frequencies and the actual temperature of the gas itself, by 
expanding the condensate cloud rapidly enough. In the condensed matter analogue, 
the typical energy of the quanta produced can in principle be even made to approach 
the relevant "Planck" scale, i.e., the point in the quasiparticle energy spectrum 
where it begins to deviate from being linear. 

Finally, on a more adventurous side, one could conceive of carrying out experi- 
ments in "experimental" cosmology, as opposed to the currently existing purely "ob- 
servational" cosmology. In such an experimental approach to matters cosmological, 
one would try to reproduce under certain specified and, in particular, well-defined 
initial conditions large-scale features of the cosmos, in the laboratory setting of 
nonstationary, inhomogeneous superfiuid gases. 
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